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The single spin asymmetry An, for large Pt single inclusive particle production in p^p collisions, 

is considered within a generalised parton model and a transverse momentum dependent factorisation 

scheme. The focus is on the Sivers effect and the study of its potential contribution to An, based 

£C) , on a careful analysis of the Sivers functions extracted from azimuthal asymmetries in semi-inclusive 

t-H ■ deep inelastic scattering processes. It is found that such Sivers functions could explain most features 

of the An data, including some recent STAR results which show the persistence of a non zero An 

CNJ t up to surprisingly large Pt values. 

U • 

O ,■ PACS numbers: 13.88.+e, 12.38.Bx, 13.85. Ni 

<\ 

C\ ■ I. INTRODUCTION 

(N; 

i— 1 1 Among the leading-twist Transverse Momentum Dependent Partonic Distribution Functions (TMD-PDFs, often 

i-Q . shortly referred to as TMDs), the Sivers distribution [l|-[3| is most interesting and widely investigated. It describes 

*~^' the number density of unpolarised quarks q (or gluons) with intrinsic transverse momentum k± inside a transversely 

O !■ polarised proton p^, with three-momentum P and spin polarisation vector S, 

<D ■ 
& f q/pt (x,k ± ) = f q/p (x,k ± ) + ^A N f q/pt (x,k ± )S-(Pxk ± ), (1) 

,. ■ where x is the proton light-cone momentum fraction carried by the quark, f q / p (x,k±) is the unpolarised TMD 
7— I ■ (k± = \k±\) and A N f q / p ^(x,k±) is the Sivers function. P = P/\P\ and k± = k±/k± are unit vectors. Notice that 
the Sivers function is most often denoted as / 13 ?(x, Aj.) [If; this notation is related to ours by Q 



t^ 



m 



2k, 



^f ; A %/pt (x,k ± ) = -—^ f£ (x,k ± ). (2) 

O ' 

A knowledge of the Sivers distribution allows a modelling of the 3-dimensional momentum structure of the nucleon [6( 
and, possibly, an estimate of the parton orbital angular momentum |7| . 

All the available information on the Sivers function has been obtained from SIDIS data, i N — »• £hX, and the study 

KJ{ • of the azimuthal distribution of the final hadron h around the 7* direction in the 7* — N centre of mass (cm.) frame. 

This analysis is based on the TMD factorisation scheme [8Ull| , according to which the SIDIS cross section is written 

as a convolution of TMD-PDFs, Transverse Momentum Dependent Fragmentation Functions (TMD-FFs) and known 

elementary interactions. Such a scheme holds in the kinematical region defined by 

P T ~ fejL ~ Aqcd « Q , (3) 

where Pt is the magnitude of the final hadron transverse momentum. The presence of the two scales, small Pt and 
large Q, allows to identify the contribution from the unintcgrated partonic distributions (Pt — k±), while remaining 
in the region of validity of the QCD parton model. The study of the QCD evolution of the Sivers and unpolarised 
TMDs - the so-called TMD evolution - has much progressed lately [8|, m, lO - fla ] , with the first phenomenological 
applications |16l - r2ll ]. 

The extraction of the Sivers functions from SIDIS data can then be performed on a sound ground. This has 
been done for the first time in Refs. (22rl26J. exploiting HERMES [13 and COMPASS [Hi data, and resulting in a 
reasonable knowledge of the Sivers functions for u and d quarks, although in a limited range of x values, x ^ 0.3. 

Much literature has emphasised the special interest and the peculiar properties of the Sivers effect. Trying to 
understand its origin at the partonic level has related the possibility of a non zero Sivers function with final |29j or 
initial [30| state interactions, respectively in SIDIS and Drell-Yan (D-Y) processes. This, in turns, induces a process 
dependence of the effect itself. The most clearcut consequence is the prediction of an opposite sign of the Sivers 



functions when contributing to SSAs in SIDIS and D-Y processes [3l|; as polarised D-Y experiments have never been 
performed so far, such a prediction has not been tested yet. Crucial information might be available in the future 
from p^p experiments at RHIC, Fermilab or from the COMPASS hadronic run at CERN, with pions colliding on a 
polarised nucleon target. The TMD factorisation scheme, valid for SIDIS processes, holds for D-Y as well, where the 
small and large scale are respectively the total transverse momentum (qx) and the invariant mass (M) of the leptonic 
pair. 

In this paper we focus on another class of puzzling results which strongly challenge our understanding of high energy 
strong interactions, that is the SSAs, usually denoted by An, measured in p^p — >• h X inclusive reactions and defined 
as: 

da^-da^ . , , +| E h d<jP'' l P^ hx 

AN = d^Td^ wlth ^ Us ¥f h • (4) 

and where f, 4- are opposite spin orientations perpendicular to the scattering plane, in the p^p cm. frame. An differs 
from the SSAs of SIDIS and D-Y processes because in such a case there is only one large scale in the process - the 
transverse momentum Pt of the hnal observed hadron - and there is no small scale related to the intrinsic motions, 
both in the distribution and fragmentation functions, which are integrated over. The TMD factorisation scheme used 
for SIDIS and D-Y processes has not been proven in this case. 

Large values of An have been measured since a long time in many different experiments. The first ones where 
at relatively low energy J32H39J . and the common expectation was that such asymmetries would vanish at higher 
energies; however, data from RHIC at ^/s = 62 A [|oj], 200 (4IH45J or even 500 (Jo. Ii?j GeV, still show puzzling non 
zero values of An ■ 

Several approaches to understanding An, within QCD and some sort of factorisation scheme, can be found in the 
literature. All of them, directly or indirectly, are related to the Sivers function or other TMDs. 

A QCD collinear factorisation formalism at next-to-lcading-power (twist-3) has been developed and used in the 
phcnomcnological studies of An [4814561 ]. In this approach the spin effect is not embedded in a spin dependent 
TMD, but the necessary phase for generating the non-vanishing SSAs arises from the quantum interference between 
an elementary scattering amplitude with one active collinear parton and an amplitude with two active collinear 
partons. The SSAs are therefore proportional to some non-probabilistic three-parton correlation functions, which are 
convoluted with product of amplitudes, rather than cross sections. These amplitudes are process dependent, while 
the three-parton correlation functions are universal. 

However, one can show that the twist-3 three-parton correlation functions have a close connection with the k±- 
moment of the TMD-PDFs; in particular the quark-gluon correlator is related to the first fci-momcnt of the SIDIS 
Sivers function [571 ]. It has been recently pointed o ut 1581) that the quark-gluon correlation functions, as obtained from 
the Sivers functions extracted from SIDIS data [22|, [23J, indeed lead to sizeable values of An, which agree in mag nitude 
with the measured ones, but with the wrong sign (the so-called sign mismatch problem). A recent analysis 1591 of the 
spin asymmetry An for single inclusive jet production in p^p collisions collected by the AjvDY experiment [601 ] does 
not show the same sign problem; however, the measured asymmetry is very small. 

An alternative, more phenomenological approach, is based on the assumption of the validity of the TMD factorisation 
also for p^p — > hX processes [H. 13. l6ll - [66l | : it generalises the usual collinear factorisation scheme (Generalised Parton 
Model, GPM) and the single inclusive cross section is written as a convolution of TMD-PDFs, TMD-FFs and QCD 
partonic cross sections. In that it adopts the same scheme which holds for SIDIS and D-Y processes with one small and 
one large scale. In this model the spin effects are included in the TMDs, which are supposed to be process-independent. 

More recently, a third approach has been proposed [67J, (68[, which assumes the TMD factorisation as in the GPM, 
but takes into account and absorbs the initial and final state interactions, i.e. the process dependence of the Sivers 
function, in the elementary interactions. In such a scheme the cross section is a convolution of process-independent 
TMDs with process-dependent hard parts; these modified hard parts are very similar in form to those in the twist-3 
collinear approach. It turns out that this modified GPM formalism leads to results and predictions opposite to those 
of the conventional GPM J67( . 

In this paper we explore the possibility of understanding the experimental results on An in p^p — > h X processes 
with the Sivers effect and within the generalised parton model of Refs. [641 66] . The first phenomenological applications 
of the Sivers effect [6ll . l62l.l64l in hadronic interactions considered the Sivers function as a free input, not constrained 
by SIDIS data. In Ref. [69| it was shown that the use of the Sivers functions, as extracted from SIDIS data, could in 
principle explain the SSAs observed at RHIC, both in size and in sign. We further pursue this study, with a careful 
analysis of the SIDIS extracted Sivers functions, with their uncertainties, and investigate whether such functions, 
assumed to be process independent, can explain the data on An, including the most recent ones. A similar study has 
been recently completed [70[ regarding the Collins effect [7l[, with the conclusion that it cannot, alone, explain all 
the available data on An ■ 



II. SIVERS EFFECT AND A N IN THE GENERALISED PARTON MODEL FORMALISM 



The Generalised Parton Model J64H661 [70] , can be considered as a natural phenomenological extension of the usual 
collinear factorisation scheme, with the inclusion of spin and k± effects through the TMDs and the dependence of 
the elementary interactions on the parton intrinsic motions. The single spin effect, da 1 ^ da^, originates from the 
TMDs; in Ref. [72| and its correction [70J ^ was shown that the only non negligible contributions to An are given by 
the Sivers TMD-PDF and the Collins TMD-FF, 

. _ [da 1 ' - rfcr 4 -] S i VCTS + [da 1 - rfo^jcoilms (f .s 

da 1 + da^ 

The Collins contribution was studied in Ref. [701, while this paper is devoted to the Sivers effect. 

In our GPM scheme the contribution of the Sivers effect to the numerator of An, for p 1 p — > h X large Pt processes, 
is given by 

[da 1 - da% ivcl . s = V f dx « dxbdz S k ^ a d 2 k ^ b d 3 p± S{p± . p c) J{p±) 5{s + i+{l ) 
*-^ ,J 16ir z x a x b z z s 



,b,c,d 



1 



x A N f a/pr (x a ,k ±a )cos(Mf b/p (x b ,k ±b )- \M»\ 2 + \M»\ 2 + \M»\' D h/c (z, P± ) , (6) 

Z L J ab-^cd 

where A N f a / p t(x a ,k± a ) is the Sivers function for parton a, Eqs. ([TJH]), which couples to the unpolarised TMD for 

parton b, f b / p {x b ,k± b ) 1 and the unpolarised fragmentation function Df l / C (z,p±) of parton c into the final observed 

hadron h. p^ is the transverse momentum of hadron h with respect to the 3-momentum p c of its parent fragmenting 

parton. J(p±) is a kinematical factor, which at 0(p±/ Eh) equals 1. For details and a full explanation of the notations 

we refer to Ref. [65| (where p ± is denoted as k±c). 

The phase factor cos(^ a ) originates directly from the k± dependence of the Sivers distribution [S-(Pxk±), Eq. (JTJ)], 

while the Mf are the three independent hard scattering helicity amplitudes defined in Ref. [65[ , describing the lowest 

order QCD interactions. The sum of their moduli squared is proportional to the elementary unpolarised cross section 
da ab^cd that ig 



da' 



±.ab-^-cd 



lit\M°\ 2 . (7) 



dt 16ns 2 2 . 

The explicit expressions of J^ |M 4 °| 2 , which give the QCD dynamics in Eq. ([5]), can be found, for all possible elementary 
interactions, in Ref. [65|. The QCD scale is chosen as Q = Pt- 

The denominator of Eq. ([?]) or ([5]) is twice the unpolarised cross section and is given in our TMD factorisation by 
the same expression as in Eq. ©, where one simply replaces the factor A N f a / p i cos((f> a ) with 2f a / p . In Ref. [69| it 
was shown that such an expression leads to results for the unpolarised cross section in agreement with data. 

We can now use the information so far available on the Sivers functions as extracted from SIDIS data and give 
some realistic estimates for the Sivers contribution to An for several single-inclusive large Pt particle production in 
p 1 p collisions. More specifically, we will consider the Sivers effect for inclusive pion, kaon, photon and jet production 
and will see how much it can contribute to the available experimental data on An ■ The analogue of Eq. (JB]) for direct 
photon and inclusive jet production will be given below. 

A. The Sivers functions in SIDIS and p^p — > hX processes 

Let us start by considering the available information on the Sivers functions and the procedure followed to obtain 
them. The first extraction - from now on denoted as SIDIS-1 fit - was presented in Ref. [22| . where the MRST01 set 
for the unpolarised PDFs [73| and the Kretzer set for the unpolarised FFs [74j were adopted. An updated extraction 



of the Sivers functions - SIDIS-2 fit - was presented in Ref. [23|. In this case, the GRV98 set for the unpolarised 
PDFs [75j and the pion and kaon FFs by de Florian, Sassot and Stratmann [70] were considered. Notice that the use 
of different PDFs does not make any relevant difference; therefore, in the following, we will consider only the GRV98 
set. 

The main features of the parameterisations adopted in those studies are the following: the analysis of SIDIS data 
is performed at leading order, 0(k±/Q), within the proven TMD factorisation approach for SIDIS, where Q is the 
large scale in the process. A simple factoriscd form of the TMD functions was adopted, using a Gaussian shape for 
their k± dependent component. For the unpolarised parton distribution and fragmentation functions we have: 



e -ki/(kD e -pi/{ P ±) 

(8) 



f q/p (x,kx) = f g /p(x) ^rp-r D h/q (z,px) = D h/q {z) ^ 

where (fcj_) and (pj_) have been fixed by analysing the Cahn effect in unpolarised SIDIS processes, see Ref. [771 ]: 

(*i) = 0.25 GcV 2 , (pi) = 0.20 GcV 2 . (9) 



The recently introduced TMD-cvolution was not taken into account, while we considered the DGLAP QCD evolution 
of the collinear factorised part. 

fq/pf 



The Sivers functions, A N f„/ p r(x, k±), have been parameterised as follows: 



—k Ilk ) 

A N f q/pt (x, fc x ) = 2Af q s (x) f q/p (x) h{k x ) 6 ' 2 " 

n \ K x) 



(10) 



where 



Af q s (x) = N s q x«« (1 - *)* K+f^^ , (11) 



with liVfl < 1, and 



a q rfi, 
OLq Pq 



h{kx) = V2-e k -±e- k ^ M \ (12) 



With these choices, the Sivers functions automatically fulfil their proper positivity bounds for any (x,kx) values. 
For the Q 2 evolution of the Sivers function, as commented above, we consider the unpolarised DGLAP evolution of 
its collinear factor f q / p (x). Notice that in the SIDIS- 1 fit we actually exploited also a different (power-like) functional 
form for h(kx), still controlled by a single parameter, leading to almost no differences in our results. In what follows 
we will only use the functional form given in Eq. (I12[) . 

In order to reduce the number of free parameters, some additional assumptions were adopted. Concerning the 
SIDIS-1 fit, we considered only u and d quarks Sivers functions, with flavour dependent a and /3 parameters. This 
amounts to a total of 7 parameters: 

N u , N d , a u , a d , p u , Pd, M. (13) 

In the SIDIS-2 fit, since we were aiming also at explaining some large kaon SIDIS azimuthal asymmetries, we 
tentatively included also the Sivers functions for antiquarks and strange quarks, u, d, s and s. To keep the number 
of parameters under control we then assumed flavour independent a and /3 parameters for the sea quarks (a soa , /3 S ea)- 
Moreover, since the large x behaviour of the Sivers function could not, and still cannot, be constrained by SIDIS data 
(see a more detailed comment below), we also assumed a single flavour independent /3 parameter, equal for quarks 
and antiquarks. This amounts to a total of 11 free parameters: 

N u , N d , N u , N s , N s , N s , a u , a d , a sca , /3, M . (14) 

Notice that even with such a choice, our complete parameterisation of the Sivers functions, Eq. (|10[) , allows for further 
differences among parton flavours, which are contained in the usual unpolarised PDFs. 

Both fits gave good results. Nevertheless it is worth stressing the main differences in the two extractions, which 
indeed play an important role in the present study. In fact, a direct use of SIDIS-1 results in the computation of 
SSAs in p^p — > hX processes for RHIC kinematics, as presented in Ref. [69|], gave very encouraging results. Notice 
that at that time the Collins effect was believed to be suppressed uM- O n the other hand, if we use the SIDIS-2 fit 
to compute the same SSAs we would get too small An values. The reason being the different j3 values coming from 
the two fits. 

More generally, as discussed in the context of the Collins SIDIS azimuthal asymmetries for the transversity distri- 
butions [70| , a study of the statistical uncertainties of the best fit parameters clearly shows that SIDIS data are not 
presently able to constrain the large x behaviour of the quark (u, d) Sivers distributions, leaving a large uncertainty 
in the possible values of the parameter j3. This is due to the limited range of Bjorkcn x values currently explored by 
HERMES and COMPASS experiments, ig < 0.3. In this respect the large ig results expected from JLab 12 GeV 
experiments will be precious |78l I79J . 



This uncertainty plays a crucial role when one tries to study the SSAs in hadronic collisions starting from the results 
obtained from SIDIS data, because the largest pion SSAs are measured at large Fcynman x values, xp > 0.3, which 
implies x > 0.3. 

To investigate the role of the Sivers effect in explaining the large value of An in p^p collisions we should therefore 
carefully explore the large x behaviour of the Sivers functions. To this aim we follow the same strategy we have 
devised in the recent study of the contribution of the Collins effect to An, the so-called "scan procedure" [70|. Here 
we summarise schematically its main steps and motivations. 

• The /3 q parameters, which control the large x behaviour of the TMDs and are largely undetermined by SIDIS 
data, play instead an important role in the computation of An, which is sizeable mainly in the large xf region. 
We can notice this explicitly by comparing, as commented above, the different implications on An of the SIDIS- 
1 and the SIDIS-2 fits. In our choice of the independent parameters it is then natural to allow for a flavour 
dependence of f3, limited, because of the relevance of the large x region, to the valence quark contributions. More 
explicitly, we only use the PDFs for u and d valence quarks in the Sivers functions (|10[) and the contribution of 
sea quarks and gluons is neglected in the sum over partons in Eq. ([()]). 

• We start the scan procedure by performing a preliminary 7-parameter [those of Eq. (fT3|) ] "reference fit" to SIDIS 
data. This reference best fit will have a total \ 2 = Xo- We then let the two parameters j3 u and (3d vary, choosing 
them in the range 0.0 — 4.0 by discrete steps of 0.5, and for each of the 81 pairs of fixed /3s we perform a new 
5-parameter fit to SIDIS data. 

• As a next step we select only those fits leading to a \ 2 such that \ 2 < Xo + ^X 2 ■ Notice that, since the reference 
fit and the scan fits have a different number of free parameters, the selection criterium is applied to the total \ 2 
rather than to the \ 2 P er degree of freedom, x\ Q v The chosen value of A% 2 is the same as that used to generate 
the error band, following the procedure described in the Appendix A of Rcf. [2j]. We find (for 217 data points) 
Xo = 270.51 and A\ 2 = 14.34. As expected from the arguments given above, all 81 fits lead to acceptable \ 2 
values for SIDIS data; this further confirms the observation that the SIDIS data arc not sensitive to the large x 
behaviour of the Sivers function. 

• We then compute, for each of the 81 selected sets, the contribution of the Sivers effect to An, according to 
Eqs. (H])-©. We do that for pion and kaon production in the kinematical regions of the STAR and BRAHMS 
experiments at RHIC. The corresponding results span the shaded areas (scan band) which are shown in the 
figures of our results. When compared with the experimental available data, the scan bands show the potentiality 
of the Sivers effect alone to account for the measured values of An in p^p — > h X processes, while preserving a 
fair description (quantified by A% 2 ) of the SIDIS data on the Sivers azimuthal asymmetry. 

• We have considered in our scan procedure all available SIDIS data, with the exception of the recent ones by the 
COMPASS Collaboration off a transversely polarised proton target [8(j. As shown in Refs. [lg, [l7| the analysis 
of these data, reaching higher Q 2 values, requires a careful use of the proper TMD evolution, which is ignored 
here, as a correct implementation of the TMD evolution in p^p — > h X large Pt processes is so far unknown. 
We have checked that the xlof °f our f^ s would be approximately (30-40)% worse for the SIDIS data including 
the proton COMPASS results and no TMD evolution. 

• We study the contribution of the Sivers effect to the SSA Am at RHIC energies only, although it might contribute 
also to the (larger) SSAs measured at lower energies [32T - [39f . The reason is that we consider only the processes 
for which our GPM and TMD factorisation can reasonably well reproduce the unpolarised cross section |69| ■ 

1. Results from the scan procedure 

Some of our results for RHIC experiments are shown in Figs. [TJ31 We have computed An adopting, as explained 
above, a single set of collinear parton distributions [75| and two different sets for the pion and kaon collinear FFs |74l . 
I76J ; the results shown correspond to the Kretzer set. Other results not shown are very similar and would not add any 
significant information. 

Let us start by considering the case of inclusive pion production. This will also help a direct comparison with the 
corresponding study on the potential role of the Collins contribution to the same observable [70[ • 

In Fig. Q] the scan band for An, as a function of Xf at fixed scattering angles, is shown for charged pions and 
BRAHMS kinematics, while in Fig. [2] the same result is given, at fixed pseudo-rapidity values, for neutral pions and 
STAR kinematics. We also give the scan band, as a function of Pt at several fixed xp values, for STAR kinematics 
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FIG. 1: Scan band (i.e. the envelope of the 81 curves obtained with the scanning procedure) for the Sivers contribution to the 
charged pion single spin asymmetries An, at ^/s = 200 GeV, as a function of if at two different scattering angles, compared 
with the corresponding BRAHMS experimental data [43J]. The shaded scan band is generated, adopting the GRV98 set of 
collinear PDFs and the Kretzer FFs, following the procedure explained in the text. 
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FIG. 2: Scan band (i.e. the envelope of the 81 curves obtained with the scanning procedure) for the Sivers contribution to 
the neutral pion single spin asymmetry An, at y/s = 200 GeV, as a function of xf at two different pseudo-rapidity values, 
compared with the corresponding STAR experimental data [4J1- The shaded scan band is generated, adopting the GRV98 set 
of collinear PDFs and the Kretzer FFs, following the procedure explained in the text. 



in Fig. [3] All these results are given at s/s = 200 GeV. We then consider the latest and interesting preliminary data 
obtained by STAR at large Pp and ^/s = 500 GeV j46(, and show our scan band in Fig. [4] for different values of xp. 

From these results we can conclude that the Sivers effect alone might in principle be able to explain the BRAHMS 
charged pion results on An in the full kincmatical range so far explored, as well as almost the full amount of STAR 
7T° data on An. This is to be contrasted with the analogous study of the Collins effect [70] , with the conclusion that 
such effect alone cannot explain the observed values of An in the medium-large xp region. 

This can be understood as follows. In the case of SSAs for neutral pion production the Collins effect suffers from 
two possible cancellations: the opposite sign between the u and d quark transversity distributions and the opposite 
sign between the favoured and disfavoured Collins FFs (necessary to build the Collins FF for 7r°); instead, for the 
Sivers effect only a cancellation between u and d flavours in the distribution sector may play a role, as it couples to 
the unpolarised TMD FF. 

The results obtained with a different choice of the fragmentation functions (the DSS set) are qualitatively very 
similar in the large xp regions. They are instead smaller in size at smaller xp, due to the large gluon contribution in 
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FIG. 3: The same as in Fig. [2j but with the STAR data plotted vs. the pion transverse momentum, Pt, for different bins in 
x F , {xf) = 0.28, 0.37, 0.43 and 0.50. 
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FIG. 4: Scan band (i.e. the envelope of the 81 curves obtained with the scanning procedure) for the Sivers contribution to 
the neutral pion single spin asymmetry An, as a function of Pt for different xf values at -^/s = 500 GeV, compared with the 
corresponding STAR preliminary experimental data at (xf) — 0.20,0.28 [4q| . The shaded scan band is generated, adopting 
the GRV98 set of collinear PDFs and the Kretzer FF set, following the procedure explained in the text. 
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FIG. 5: Scan band (i.e. the envelope of the 81 curves obtained with the scanning procedure) for the Sivers contribution 
to the kaon single spin asymmetry An, as a function of xf, at y/s = 200 GeV and a fixed scattering angle, compared with 
the corresponding BRAHMS experimental data [43]. The shaded scan band is generated, adopting the GRV98 set of collinear 
PDFs and the Kretzer FF set, following the procedure explained in the text. 

the leading order (LO) DSS fragmentation functions. They are not shown here. 

For completeness we consider the case of SSAs for kaon production as measured by the BRAHMS collaboration. 
Our results for the scan band, compared with the available data, are shown in Fig. [5] 

The results in Figs. [T][5] show that the Sivers effect alone, as computed in our GPM scheme, Eqs. ((I])-©, and based 
on the Sivers functions extracted from SIDIS data and assumed to be universal, can be large enough to explain alone 
the pion SSAs An observed at RHIC. One should not forget that, indeed, the phenomenology of the Sivers effect was 
originally generated in the attempt to explain the large values of A n observed by the E704 Collaboration [l|, y, [fLU] . 
However, the amount of uncertainty in the scan bands, due to lack of precise SIDIS data at large x, is still much too 
large to draw any definite conclusions. 

A full understanding of the SSAs in inclusive p^p — > hX processes should also take into account the contribution of 
the Collins effect, which might be small, but not entirely negligible. Rather than addressing the issue of a best fit of 
SIDIS + An data with Collins and Sivers effects, which is premature at this stage, we now adopt a more pragmatic 
attitude. We wonder whether, among the 81 sets of parameters which build up the possible results on An contained 
in the scan bands, we can find some which give a good description of all the data. 



2. Results with a selected set of parameters and its statistical uncertainty bands 



Among the full set of curves produced by 
of An (actually one could find more than 
corresponding statistical error band. Our 
■k^ data vs. xp at fixed angles, for STAR 
at different xp values, and for BRAHMS 
corresponding values of the parameters are 
a set of Sivers functions for u and d quarks 



the scan procedure, we have isolated the set leading to the best description 
a single set); we have then evaluated, as in Appendix A of Ref. [23|], the 
results are presented in Figs. El El El an d El respectively for BRAHMS 
77° results vs. xp at fixed pseudo-rapidities, for STAR tt° results vs. Pp 
K^ data vs. xp at a fixed angle, all of them at y/s = 200 GeV. The 
given in Table 1. From these results one can see that it is possible to find 
which, while describing well the SIDIS data, can also describe fairly well, 



TABLE I: Our chosen set of 7 parameters, Eq. (|13|l . fixing the u and d quark Sivers distribution functions, according to 
Eqs. (|10I12I) . Among the 81 sets of the scan procedure, this set gives the best description of the Am data. The corresponding 
total value of \ 2 for the 217 SIDIS data points is 273.2, which is very close to the best value Xo — 270.5 of the reference set. 
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FIG. 6: The Sivers contribution to the charged pion single spin asymmetry An, compared with the corresponding BRAHMS 
experimental data at two fixed scattering angles and y/s = 200 GeV J4J|. The central lines are obtained adopting the GRV98 
set of collinear PDFs and the Kretzer FFs, with the Sivers functions as in Eqs. (|10|| - (|12|l with the parameters given in Table 1. 
The shaded statistical error bands are generated applying the error estimate procedure described in Appendix A of Ref. [23(1 . 
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FIG. 7: The Sivers contribution to the neutral pion single spin asymmetry An, compared with the corresponding STAR 
experimental data at two fixed pion rapidities and y/s = 200 GeV [4J] . The central lines are obtained adopting the GRV98 set 
of collinear PDFs and the Kretzer FFs, with the Sivers functions as in Eqs. (I10|l (|12|) with the parameters given in Table 1. 
The shaded statistical error bands are generated applying the error estimate procedure described in Appendix A of Ref. [23|] . 



alone, the SSAs for pion production, as measured both by BRAHMS and STAR Collaborations at 200 GeV. 

The preliminary STAR data at 500 GeV [46[ deserve a dedicated comment. Quite surprisingly, they show values of 
An of the order of few percents, with a flat behaviour as a function of Pt at fixed xp, up to Pt — 7 GeV. Such a trend 
is well reproduced by our set of chosen best parameters; however the computed magnitude of An is smaller than data, 
as shown in Fig. 1101 left plots. Being the asymmetry so small, we have also computed the Collins contribution to An, 
following Ref. [701 ]. It turns out that, for some sets of the parameters, the Collins contribution has a similar trend and 
magnitude as the Sivers one, as shown in Fig. 1 101 right plots. Then, an appropriate sum of the two contributions, 
according to Eq. ([5]), might well explain also this new puzzling data. 

Another cautious comment about the STAR data on An at 500 GeV concerns the large value of their QCD scale, 
Q 2 — P T . As we noticed for the COMPASS proton data, at such values the TMD evolution might play an important 
role. Our results should then be taken as an indication in favour of a combined Collins + Sivers effect, rather than a 
proof. Qualitatively, one expects from TMD evolution an increase of the average (fcj_) value of the Sivers distribution, 
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FIG. 8: The same as in Fig. [7J but with the STAR data plotted vs. the pion transverse momentum, Pt, for different bins in 
x F , {x F } = 0.28, 0.37, 0.43 and 0.50. 



which would help increasing the corresponding value of A 



N- 



B. SSAs for p^p — > jetX and p^p -47I processes 

In these processes no fragmentation mechanism is required, so that, within the GPM and the TMD factorisation 
approach, one can access directly the spin and k± properties of the partonic distributions. After integration over the 
intrinsic azimuthal phases only the Sivcrs effect survives, which is then best studied in these processes, as discussed, 
e.g., in Refs. [64|, Hl(. Notice that, for the same reasons, the SSAs for inclusive jet or photon production can be used 
to test the process dependence of the Sivers functions in a modified generalised parton model with inclusion of initial 
and final state interactions [671 182| or within the twist-3 approach [591 ] . 

The numerator of An for the inclusive jet production can be obtained from Eq. ^ simply replacing the TMD 
fragmentation function, D !l / C (z,p±), with a factor 6(z — 1) S 2 (p±) (and identifying now the final hadron momentum, 
p h , with the jet momentum p c = Pj ct ). More explicitly the numerator of An for inclusive jet production reads 



[d^-da^ c T tX ■ 



a,b,c,d 



dx a dxb 

16 7T 2 X a Xb s 



d 2 kj_ a rf 2 fcj_ 6 S(s + i+u) 



x A n f a/p r(x ai k±a) cos(0 a ) f b/p (x b , k ±b ) - \\M^ + \M»\ 2 + |M 3 U | 



(J 1 2 



ab—tcd 



(15) 



0.1 



0.05 



- 



-0.05 



— i 1 1 1 1 1 1 1 1 1- 

K + — — ' 



9 = 2.3 



_i i i i i i i i i_ 



stat. bands 



n 



— i 1 1 1 1 1 1 1 1 1— 

K — — ■ 



*1^ 



-j i i i i i i i i_ 



0.2 0.3 0.4 0.2 0.3 0.4 



*F *F 

FIG. 9: The Sivers contribution to the charged kaon single spin asymmetry An, compared with the corresponding BRAHMS 
experimental data at a fixed scattering angle and *Js = 200 GeV [43] ■ The central lines are obtained adopting the GRV98 set 
of collinear PDFs and the Kretzer FFs, with the Sivers functions as in Eqs. (1101) (I12[) with the parameters given in Table 1. 
The shaded statistical error bands are generated applying the error estimate procedure described in Appendix A of Ref. [23] . 
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FIG. 10: Left panels: the Sivers contribution to the it single spin asymmetry An vs. the pion transverse momentum Pt, for 
different bins in xf, compared with the corresponding STAR preliminary data at y/s = 500 GeV and (xf) = 0.20, 0.28 [46[ [. 
The central lines are obtained adopting the GRV98 set of collinear PDFs and the Kretzer FFs, with the Sivers functions as 
in Eqs. (|10[) - ()12[) with the parameters given in Table 1. The shaded statistical error bands are generated applying the error 
estimate procedure described in Appendix A of Ref. [2a |. Right panels: the Collins contribution to the same An, computed 
according to Ref. [7Q], choosing the Collins functions, among those of the scan band, which give the maximum contribution. 
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functions as in Eqs. (|10|1 (|12p with the parameters given in Table 1. The shaded statistical error band is generated applying 
the error estimate procedure described in Appendix A of Ref. [23)]. Right panel: the same estimate as in the left panel for a 
direct photon, rather than a jet, production at y/s = 200 GeV and n = 3.5. 
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FIG. 12: Our computation of the unpolarised cross section for jet production vs. the jet energy, at y/s = 510 GeV and fixed 
pseudo-rapidity r\ = 3.25, compared with AjvDY data [47| . 



Notice that the elementary hard scattering interactions are exactly the same as those for the inclusive hadron pro- 
duction and the jet, at LO, is identified with the final parton c. 

Concerning the direct photon production the basic partonic processes are the Compton process gq{q) —> jq (q) 
and the annihilation process qq — > jg. In this case one can formally use the above equation replacing the partonic 
unpolarised cross section, Eq. (|7|), with the corresponding one for the process ab — > ^ d (see also Ref. |64|). 

No SSA data are so far available for direct photon production, while very recently some preliminary data for 
inclusive jet production have been released by the A/yDY Collaboration at yfs = 500 GeV [43]. The values measured 
for An are very tiny, but very precise and might indicate a non zero asymmetry. 

In the left plot of Fig. [TT]we show our estimate, based on the chosen best set parameters of Table 1, for Aj\t(xf) in 
P^P ~^ jet-X" processes at a fixed pseudo-rapidity value and y^s = 500 GeV, and compare it with the A^DY data [471 ]. 
In the right plot we give our corresponding estimates for An(xf) in p^p -> 7X processes at a fixed pseudo-rapidity 
value and yfl = 200 GeV. 

For consistency, in Fig. [T^] we compare our (leading order) computation of the cross section for jet production as 
given by Eq. ([T5j) where we replace the factor A N f a / p t cos(0 a ) with f a /p, with the A^DY data at y/s = 510 GeV and 
fixed pseudo-rapidity r\ = 3.25. 
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III. CONCLUSIONS 

We have explored, within our generalised parton model scheme which assumes TMD factorisation and the uni- 
versality of the parton distribution and fragmentation functions, the contribution of the Sivcrs effect to the single 
spin asymmetry An as measured by RHIC Collaboration experiments. The Sivers functions are the same as those 
extracted from SIDIS processes. A similar analysis was performed in Ref. [7(|, concerning the Collins effect. 

Our results, limited to the contribution of the valence quarks, to the SIDIS cases where the TMD evolution is not 
expected to be relevant and to the pp — > tt X large Pt processes for which we can well reproduce the unpolariscd 
cross sections, are rather encouraging. For most pion data the Sivers effect alone could explain the observed values of 
An in magnitude and, in particular, in sign. This is in contrast to other approaches, also related to the Sivcrs effect, 
which seem to have severe problems [58l. l67| in explaining the sign of the observed An- 

We have performed our analysis by varying the parameters of the Sivers functions which are not well fixed by the 
SIDIS data, due to their limited kinematical range, obtaining the so-called scan bands. Then, among the explored 
sets of parameters, wc have shown that a particular choice gives a reasonable description of the An data. We have 
also used such a set to compute estimates for SSA in p'p — > jet X and p^p — > 7 X processes. Such measurements will 
further allow to discriminate between our approach and others. 

While encouraged by the results of our analysis we should avoid making definite conclusions at this stage. This 
work shows that the GPM TMD factorisation scheme could explain at the same time the main features of the SSAs 
measured in SIDIS and hadronic processes. While such a scheme is well justified for SIDIS it can, so far, only be 
considered as a phenomenological model for hadronic processes, which needs further confirmation or disproval from 
data and further theoretical work. Our choice of the sets of parameters given in Table 1 is not meant to be interpreted 
as the final best set of Sivers functions. A full analysis, including TMD evolution, of all data involving the Sivcrs 
effect - i.e. all the SSAs in several different processes - would require much more attention and work. 
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